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In this work, we present a work in progress towards an efficient and economical computational 
module which interfaces between Cauchy and characteristic evolution codes. Our goal is to provide 
a standardized waveform extraction tool for the numerical relativity community which will allow 
CCE to be readily applied to a generic Cauchy code. The tool provides a means of unambiguous 
comparison between the waveforms generated by evolution codes based upon different formulations 
, of the Einstein equations and different numerical approximation. 
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PACS numbers: 04.20Ex, 04.25Dm, 04.25Nx, 04.70Bw 



INTRODUCTION 



The study of gravitational waves will brighten unexplored features of our universe that are otherwise invisible to 
conventional astronomy and will increase our knowledge about the very nature of time and space Gravitational 
wave detectors are already operating, and results from the first LIGO and Virgo collaboration were recently published 
in Nature 0. The signal predicted by numerical relativity will provide a template bank used for filtering the noise, 
indispensable to the success of gravitational wave detectors such as LIGO, Virgo, and LISA. The current sensitivity 
levels of the detectors will be improved substantially in next-generation detection estimated by 2015. Although 
existing simulations are sufficiently accurate for populating the parameter space in current searches of ground-based 
detectors, the new generation of advanced detectors will be 10-15 times more sensitive by 2015. 

Ideally, the emitted gravitational wave signature should be extracted at spatial or null infinity. However, most 
present codes impose artificial, finite outer boundaries and are performing the waveform extraction at finite radius. 
■ This method introduces systematic errors, especially for higher modes, which is the main obstacle in reaching the 
desired accuracy. With the exception of Pretorius |3||, who uses coordinates that compactify spatial infinity, all the 
other codes use a computational domain with a finite outer boundary and Sommerfeld-like approximate outer boundary 
conditions must be imposed, which introduce errors in the computation of gravitational waves. The choice of proper 
boundary conditions is complicated by gauge freedom and constraint preservation Q. The emitted gravitational 
wave signature is calculated at finite distance, using either the Newman-Penrose Weyl scalar ip 4 or the odd and 
even parity functions Q + , Q x in the Zerilli-Moncrief formalism Q. The strain h of the wave used in detection is 
obtained performing one time integration from the Zerilli-Moncrief multipoles, and two time integrations from the 
Newman-Penrose curvature. The waveform is affected by gauge ambiguities which are magnified by the integration Q. 

Cauchy-characteristic extraction (CCE) [8j, which is one of the pieces of the CCM strategy 9], offers a means 
to avoid the error introduced by extraction at a finite world-tube. In CCE, the inner world-tube data supplied by 
the Cauchy evolution is used as boundary data for a characteristic evolution to future null infinity I + , where the 
waveform can be unambiguously computed by geometric methods. This characteristic initial-boundary value problem 
based upon a timelike world-tube [l(| has been implemented as a mature evolution code, the PITT null code [TTI, ll^. 
which incorporates a Penrose compactification of the space-time. By itself, CCE does not use the characteristic 
evolution to inject outer boundary data for the Cauchy evolution, which can be a source of instability in full CCM. 

The PITT code has been tested to be second order convergent in a wide range of testbeds extending from the 
perturbative regime [l3j to highly nonlinear single black hole spacetimes [l2j]. However, in cases which require high 
resolution, such as the inspiral of matter into a black hole, the error in CCE has been a troublesome factor in the 
postprocessing phase [14j . This has motivated a recent project [l5[ to increase the accuracy of the PITT code. Other 
results achieved with previous versions of the PITT have been recently reported [l6l [r7| . Recently, the code underwent 



major improvements and corrections to previous versions, to improve accuracy and convergence 181 ]. 

Here we test this improved version of CCE on a realistic application involving a Cauchy evolution of the inspiral 
and merger of two equal mass non-spinning black holes. We use the same code specifications described in [l5j except 
that the accuracy of angular derivatives has been increased to a 4th order finite difference approximation. The results 
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presented here are a work in progress towards our goal to develop CCE as a reliable and accurate waveform extraction 
tool for the numerical relativity community. This paper addresses the first two objectives: 

• To create a robust and flexible interface between a binary black-hole Cauchy evolution code and a characteristic 
code for wave extraction at infinity. 

• To prove the robustness of the interface by performing precise computations of gravitational waveforms at 
infinity from binary black-hole, using this Cauchy-characteristic extraction approach. 

We construct an interface that takes the Cartesian data from a Cauchy evolution and converts it into boundary 
data on a spherical grid for the characteristic evolution. The data are evolved to future null infinity, where it is 
used to compute the gravitational waveform. The flexibility of the interface is due to implementation of a spectral 
decomposition of data. This implementation has been tested with a realistic application involving a binary black-hole 
inspiral. In Sec.|TI]we review the formalism underlying CCE, including enough details of the patching, evolution and 
extraction, to make clear the difficulties underlying the calculation of an accurate waveform at I + . In Sec. IIIII we 
briefly describe the initial data for Cauchy and characteristic evolution. In Sec. lIVi we present the details of the CCE 
interface which allows the data from a Cauchy evolution to be used as boundary data on an inner worltube for a 
characteristic evolution to I + , where the waveform is extracted. In Sec. El we test the CCE interface by extracting the 
waveform from a Cauchy evolution of a binary black-hole inspiral and merger, and by comparing it to the waveform 
obtained by other standard method in current practice. 



Characteristic data are provided by the Cauchy evolution on a world-tube WT, free initial data being given on the 
initial null hypcrsurface Mx, which sets the metric on the entire initial cone (fig. 1). 



FIG. 1: Cauchy and characteristic evolution are patched in the vicinity of a world-tube WT, embedded in Cauchy evolution 

The metric data from a Cauchy evolution are interpolated onto a timelike inner world-tube to extract the boundary 
data for the characteristic evolution. The characteristic evolution is embedded into the Cauchy evolution and is 
extending to future null infinity T + , where the waveform can be unambiguously computed using the geometric methods 
developed by Bondi et al , Sachs and Penrose [2l| . The extraction process involves carrying out the complicated 
Jacobian transformation between the Cartesian coordinates used in the Cauchy evolution and the spherical null 
coordinates used in the characteristic evolution (the full details are given in [22j].) 



II. THE CCE FORMALISM 



A. Cauchy-Characteristic Patching 





B. Characteristic Evolution 



The characteristic formalism is based upon a family of outgoing null hypersurfaces, emanating from some inner 
world-tube, which extend to infinity where they foliate I + into spherical slices. We let u label these hypersurfaces, 
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x A (A = 2, 3) be angular coordinates which label the null rays and r be a surface area coordinate, (fig. 2). 



Timelike 
wurklmbe 




(9.0) 



FIG. 2: Ongoing null hypersurfaces emanating from the world-tube and extending to X + 
In the resulting x a — (u,r,x A ) coordinates, the metric takes the Bondi-Sachs form (l9l. l20l| 

ds 2 = -^e 2 ^-r 2 h AB U A U B ^jdu 2 -2e 2p dudr 
- 2r 2 h AB U B dudx A + r 2 h AB dx A dx B , 



(2.1) 



where h A B is the Bondi-Sachs conformal 2- metric with h AB hsc = &c- ^he c °de introduces an auxiliary unit 
sphere metric q A B, with associated complex dyad q A satisfying q A B = \ {^aQb + QaQb)- For a general Bondi-Sachs 
metric, the full nonlinear h A s is uniquely determined by the dyad component J = h A sq A q B /2, since the other dyad 
component K — h A sq A q B '/2 is constrained by the determinant condition 1 = K 2 — J J. The spherically symmetric 
case characterized by J = 0. We introduce the spin-weighted fields U = U A q A and Q = Q A q A , where 

Q A = r 2 e- 2 ^h AB U B r . (2.2) 

[23l | for further details regarding numerical 



as well as the (complex differential) operators 9 and 9. Refer to 
implementation. The auxiliary variables 

v = 3J, B = Qfi, k = dK 



(2.3) 



are also introduced to eliminate all second angular derivatives. In certain applications this has been found to give rise 
to increased accuracy by suppressing short wavelength error [2~ij . 

In this formalism, the Einstein equations = decompose into hypersurf ace equations, evolution equations and 

l26j . the hypersurface equations 

(2.4) 
(2.5) 
(2.6) 



conservation conditions on the inner world-tube. As described in more detail in 12 1 
take the form 

= Np, 
U r = r- 2 e 2p Q + Nu, 
(r 2 Q). r = -r 2 (dJ + m) ir + 2r 4 d(r- 2 p) r + N Q 
1 o«_ «„= a 1 



where 123 



n = 2K -mx + -{B 2 J 



+ -r- 2 (r 4 (m + W)) r + N v , 



9J) + — mJQj-Bjdj) 

AK 



(2.7) 



(2- 



is the curvature scalar of the 2-metric h A s- Those equations have a hierarchical structure in [J, ft, Q, U, V] such 
that the right hand sides, e..g. Np[J] only depend upon previous variables and their derivatives intrinsic to the 
hypersurface. 
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The evolution equation takes the form 

2 (rJ) ur - (r^V (rJ) r ) r = -r- 1 (r 2 W) r + 2r~ 1 e^ 2 e p - (r^V) r J + Nj, (2.9) 

where, Np, Nu, Nq, Ny and Nj are nonlinear terms which vanish for spherical symmetry. Expressions for these 
terms as complex spin- weighted fields and a discussion of the conservation conditions are given in [8|. 

The characteristic Einstein equations are evolved in a domain between an inner radial boundary at the interior 
world-tube, and an outer boundary at future null infinity. The characteristic evolution code implements this formalism 
as an explicit finite difference scheme, based upon the compactified radial coordinate 

£ = -^— (2-10) 

R E + r 

so that £ = 1 at I + . Here Re is a parameter based upon the extraction world-tube, which in the CCE module is 
chosen as the radius of the extraction world-tube, as determined by R 2 = 8ijX l x^ in terms of the Cartesian coordiates 
x % used in the Cauchy evolution code. The boundary data for J, (3, U, Q, and V on the world-tube supply the 
integration constants for a radial numerical integration of the hypersurface Einstein equations. The finite difference 
scheme for integrating the hypersurface and evolution equations is based on the marching equation for a spherically 
symmetric scalar field $: 

1 f fV\ $ 

$n - $w - *e + $s = -- J y—j -dudr (2.11) 

where the point N is the "new" point in the evolution scheme, and V is defined by the spherically symmetric version 
of the Bondi-Sachs metric given above. The evolution scheme in the full gravitational case used to determine the 
metric at the next point on the null hypersurfaces is modeled after this example (see [l2|, [24j for details) . 




FIG. 3: The null parallelogram WSEN used to determine the field values at point N, as described by (|2.11[) . 

C. Gravitational Radiation Calculation 

The theoretical derivation of the waveform at infinity is carried out in terms of an inverse surface-area coordinate 
1= 1/r, where £ — at I + . In the resulting ir/ 1 = (u 7 £ 7 x A ) coordinates, the physical space-time metric (|2.ip has 
the conformal compactification g^ u = £ 2 g^ u , where g^ v is smooth at I + and takes the form [To| 

g^dx^dx" = - (e 2f! V£ 3 - h AB U A U B ) du 2 + 2e 2f3 dud£ - 2h AB U B dudx A + h AB dx A dx B . (2.12) 

As described in [l5j], the Bondi news function N(u, x A ) and the Newman-Penrose Weyl tensor component ^{u, x A ) = 
linv^oo 7-04 which describe the waveform are both determined by the asymptotic limit at I + of the tensor field 

£»u = J(V M V„ - -^V Q V Q )£ (2.13) 
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constructed from the leading coefficients in an expansion of the metric in powers of £ 



h AB = H AB +£c AB + 0(£ 2 ), (2.14) 

/3 = H + 0(£ 2 ), (2.15) 

U A = L A + 2£e 2H H AB D B H + 0(£ 2 ), (2.16) 

f 2 V = D A L A +£(e 2H K/2 + D A D A e 2H ) + 0{£ 2 ), (2.17) 



where 1Z and D A are the 2-dimensional curvature scalar and covariant derivative associated with H AB . 

The expansion coefficients iJ, H ABl c AB and L A (all functions of u and x A ) completely determine the radiation 
field. Before the gravitational radiation is calculated from the metric in the neighborhood of Z + , it is necessary to 
determine the conformal factor uj relating H AB to a unit sphere metric Q ABl i.e. to an inertial conformal Bondi 
frame [l(| satisfying 

Q AB =u 2 H AB . (2.18) 

The news function N(u,x A ) is directly computed by the code in terms of the computational coordinates (u,x A ), as 
opposed to the inertial coordinates (u,y A ) on T + corresponding to an idealized distant observatory. The transforma- 
tion to inertial coordinates proceeds first by introducing the conformally rescaled metric g^ v = uj 2 g^ v in which the 
cross-sections of I + have unit sphere geometry, in accord with (12.18[) . Then the rescaled null vector ri/' = u; -1 /^ 
is the generator of time translations on Z + , i.e. n^du — d^. The inertial coordinates thus satisfy the propagation 
equations 

h»d„u = LJ, h^d^y A =0 : (2.19) 

where n^d^ — e~ 2H (d u + L a 8 x a) in terms of the computational coordinates. The inertial coordinates are obtained 
by integrating (I2.19[) . thus establishing a second pair of stereographic grid patches corresponding to y A . Then the 
news function is transformed into N(u,y A ). 
The Bondi news function N is given by (|2.20[) . 

N = ^e- 2l5 u- 2 e- 2H F A F B {(d u + £ L )c AB ~ \c AB D c L c + 2u J D A [u J - 2 D B (u J e 2H )}}, (2.20) 

where £l is the Lie derivative with respect to L A . The Newman-Penrose Weyl tensor component \& is given by (12 ,21[) 

* = \u- 3 e~ 2lS h»F A F B (d^AB - d A ±^ B - f« B £ Aa + T^st^ | I+ . (2.21) 

In the inertial Bondi coordinates, the expression for the news function (|2.20j) reduces to the simple form 

N = ^Q A Q B d u c AB) (2.22) 

and (|2.21l) reduces to the single term 

* = \Q A Q B d 2 u c AB = dldiJ\ I+ . (2.23) 

This is related to the expression for the news function in inertial Bondi coordinates by 

* = d u N. (2.24) 

Equation (|2.24p holds true in the linearized approximation of the Einstein equations. In the nonlinear case, the full 
expression for news and 'J must be used in the code. This introduces additional challenges to numerical accuracy due 
to high order angular derivatives of uj and large number of terms. 



III. INITIAL DATA 
A. Initial Cauchy Data 

For the Cauchy evolution we used the LazEv code [13, HH along with the Cactus framework [2§| and Carpet [3(3] 
mesh refinement driver. LazEV is an eighth-order-accurate finite-difference code based upon the Baumgarte-Shapiro- 
Shibata-Nakamura (BSSN) formulation l3ll . [32| of Einstein's equations, which deals with the internal singularities 
by the moving puncture approach (27l . |33| . Our simulation used 9 levels of refinement with finest resolution of 
h = Af/53.76, and outer Cauchy boundary at 400A/. The initial data consisted of a close quasicircular black-hole 
binary with orbital frequency Mfl = 0.050, leading to more than a complete orbit before merger (see [Uj])- We output 
the metric data on the extraction world-tube every At = M/32. 
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B. Initial Characteristic data 

The initial data for the characteristic evolution consist of the values of J on the initial hypersurface u — 0. One way 
of supressing incoming radiation in the data would be to set the Newman-Penrose Weyl tensor component Wo = 
on the initial null hypersurface. For a perturbation of the Schwarzschild metric, this condition implies no incoming 
radiation in the linearized approximation. However, in order to avoid shocks arising from incompatibility with the 
Cauchy data on the extraction world tube £ = £e (with £ given by 12.101) . we also need to require that J and d^J 
are continuous. In the linearized approximation, the condition that Wo = implies that 9|J = 0. The combination 
of those requirements leads to J = J\^ E + (c\ J)\$ E (£ — £e), which would imply that J ^ at I + . For technical 
simplicity we avoid this complication by initializing J according to 

(£- 1) 

J = J| ^|M)' (3 - 1} 

which matches the Cauchy data and the derivatives at £ = £g and is consistent with asymptotic flatness. Since this 
choice of J vanishes at infinity, the initial slice of X + has a unit sphere metric so that the conformal factor has the 
simple initialization u(0,p, q) = 1. 

Given the initial data (|3.ip . this leads to complete knowledge of the metric on the initial null cone. Then (|2.9[) gives 
an expression for J >uri which is used to determine J on the "next" null cone, so that the process can be repeated to 
yield the complete metric throughout the domain, which extends to I + . 

IV. COMPUTATIONAL INTERFACE 

We have designed an interface that takes Cartesian grid data from a Cauchy evolution and converts it into boundary 
data for characteristic evolution on a spherical grid extending to I + . We treat each component g^t, x % ) of the Cauchy 
metric as a scalar function in the x l Cartesian coordinates which are used in the 3 + 1 evolution. 

In order to make the interface as flexible as possible for future development as a community tool for waveform 
extraction, we have based it upon a spectral decomposition of the Cauchy data in the region between two world tubes 
or radii R — R% and R = R2, where R = 5ijX l x^ is the Cartesian coordinate radius. Then at a given time t = T, 
we decompose g^^T, x l ) in terms of Tchebychev polynomials of the second kind Uk(R) and spherical harmonics 
Yi m (6,4>), where (9, 4>) are related to x l /R in the standard way. The Tchebychev polynomials are conventionally 
defined as functions Uk(r) on the interval — 1 < r < 1. Here we map them to the interval Ri < R < R 2 by the 
transformation 

2R — Ri — R 2 
t(R) = . 

R 2 — R\ 

where the extraction shell thickness is determined by the number kuax of Tchebychev polynomials used. (In tests 
of binary black holes with mass M we use a relatively small range R 2 — R\ — 10M, a larger value of kuax would be 
needed if the range were expanded). Thus, for R\ < R < R 2 , we expand 

g^{T lX l ) =J2 C ^[kim]Uk(R)Yi m (e,c(>)- (4-1) 

klm 

For the applications to waveform extraction given in this paper, it is sufficient to consider I < I Max, where I Max = 
6, and k < kMax, where kuax — 6. The coefficients C^ny^ then allow us to reconstruct a spherical harmonic 
decomposition of each component of the Cauchy metric on the extraction world-tube R = Re, i.e. 

g ^[im\{T , Re) = / ] Cfj, v \ki m ]Uk(RE)- (4.2) 

This decomposition is carried out at a sequence of Cauchy time steps Tjf = Tq + NAT, where AT is chosen to be 
much smaller than the characteristic time scale of the problem but, for purposes of economy, larger than the time step 
used for the Cauchy evolution. A fifth-order polynomial interpolation is carried out locally over the Tjv to provide 
characteristic boundary data at any time t in analytic form. 

The extraction module also requires the derivatives dtg^ v and d^g^ at the extraction world-tube. The t-derivative 
is constructed by a fourth-order- accurate finite-difference stencil using the surrounding Cauchy times t = T^r. The 
i?-derivative is obtained analytically, at each time level Tpj, by differentiation of the Tchebychev polynomials. 
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The spherical harmonic interpolator from the Cartesian to the spherical coordinates is part of the extraction module, 
but its resolution is controlled by the Cauchy evolution. 

The stereographic coordinates x A = (q,p) used to label the outgoing null rays in the Bondi metric are matched 
to the spherical coordinates (9, (/>) induced by the Cartesian Cauchy coordinates on the extraction worldtube by a 
standard transformation, using the conventions in [23j . The value of the surface-area coordinate r in the Bondi- 
Sachs metric is obtained on the extraction world-tube from the 2-determinant of the Cartesian metric on the surfaces 
t = Tat, R = Re- As a result the radius of the Bondi coordinate r ^ const on the extraction world-tube. The metric 
has to be calculated at a common value of the surface coordinate r, because the original Cauchy extraction was at 
constant R. In order to make this calculation possible, the transformation from Cartesian coordinates (t, x l ) to Bondi- 
Sachs coordinates (u,r, x A ) is carried out via an intermediate Sachs coordinate system (u,X,x a ) [20j where A is an 
affinc parameter along the outgoing null rays. The afhne freedom allows us to set A = on the extraction world-tube 
R = Re- After carrying out the Jacobian transformation from (i, x l ) to (u, A, x A ), the Cartesian metric and its first 
derivatives at the extraction world-tube provide a first-order Taylor expansion in A (about A = 0) of the null metric 
in Sachs coordinates. The corresponding Taylor expansion of the metric in Bondi-Sachs coordinates then follows from 
the computed value of r and d\r at A = 0, which are obtained from the 2-determinant of the Cartesian metric. In 
order to obtain a first-order Taylor expansion for the Bondi metric variable f3, the hypersurface equation (|2.4[) must 
be used to evaluate d r /3 at the extraction world-tube. All other metric variables are then initialized consistent with 
second order accuracy. Taylor expansions are also needed to start up the radial integration equations for the auxiliary 
variables (|2.3[) used to convert angular derivatives to first-order form. These expansions are obtained from applying 
the 3-operator to the Taylor expansion of the underlying metric. This is a complicated process because the 3 operator 
intrinsic to the A = extraction world-tube is not the same as the 3 operator intrinsic to the r = const Bondi spheres 
(see [lH for a discussion). The low order intermediate Taylor expansions limit the accuracy of the result. A new 
approach that avoids entirely the use of the Taylor expansion and gives better accuracy is presented in [l8j . In the 
original approach, used for this paper, the resulting Taylor expansion of the evolution variables is used to fill the 
points of the Bond-Sachs grid to start the integration of the characteristic hypersurface and evolution equations (|2.4[) 
- (|2.9I) . The integration proceeds from the extraction world tube to I + on a radial grid based upon the compactificd 
cc-coordinate (|2. 10|) . 

Domain of dependence considerations place a constraint between the characteristic time step Au and the size of the 
characteristic grid analogous to the CFL condition for the Cauchy evolution. For an estimate, consider the Minkowski 
space case with the conformally rescaled metric 

ds 2 = — — „^ du 2 — —dud^ + qAsdx A dx B (4-3) 



'E 

where the unit sphere metric takes the form 



R 2 F Re 



q AB dx A dx B = ———(dp 2 + dq 2 ). (4.4) 
1 + p + q 



The past light cone is determined by 



du_ _ -<£ - y/di 2 + (1 - t) 2 q AB dx A dx B 



Re (1 - 

For typical characteristic grid parameters, Ap = Aq = A£/4, the resulting restriction is 

\Au 



Re 



< 8A£ (4.6) 



For a Cauchy simulation of a binary black-hole system of total mass M with timestep At — M/32 (sufficient to 
describe the typical frequencies of a binary system), (|4.6I) leads to 

M 

~ < A£, (4.7) 



256R E 

for the choice of characteristic timestep Am = At. The corresponding number of radial gridpoints must roughly satisfy 
N% < 128Re/M. This places no limit of practical concern on the resolution of the characteristic evolution even for 
the small extraction radius Re = 20M. Thus, for purposes of CCE, there are no demanding CFL restrictions. 

The interface was debugged and calibrated using the analytic Schwarzschild metric in Kerr-Schild coordinates (t, x z ), 

2m 

g^u = r\p,v -\ fc M fc„, (4.8) 



where — (— l,x l /r). 
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V. RESULTS 

We present results for the characteristic extracted waveform either in terms of ^, related to the Bondi news by 
VP = d u N in the linearized regime, or, when comparing to the perturbative waveform, in terms of the Newman-Penrose 
component ip 4 . The relationship between the Cauchy and the characteristic waveforms is: (R — 2M)ip 4 = —2^>. We 
decompose the signal in I = 10 spherical harmonic modes but, for illustrative purposes, we concentrate on the dominant 
(2,2) and sub-dominant (4,4) modes. The Cauchy data were given at the extraction radii R — 20M, 50M, 100M. 
The relationship between the Cauchy radius R and the characteristic world-tube radius Re is: Re/ R = 1 + 1/R + 
l/(4i? 2 ). The characteristic extraction module was run with the following specifications: angular gridpoints = radial 
gridpoints = 60, 120, 240, and timestep Au = 8At, 4At, 2At, where At = M/32. The test was run until t/M = 385, 
using 4*' t -order accurate angular derivatives, on stereographic patches with circular boundaries and angular dissipation 
ejx — 0.001 (see [l5| for details on how the angular dissipation is added to the evolution equation 12.9ft . The results 
are shown for the highest resolution. Table U gives the convergence rates for the world-tube variables obtained with a 
small extraction radius Re = 20M at a time corresponding to the peak of the signal (t w 200M). The rates are given 
for the real and the imaginary part. All quantities are very close to second order convergent, including J x , which is 
the term which determines the waveform. 

TABLE I: Convergence rates of the I — 2, m = 2 mode for the metric variables measured near the peak of the signal (t ~ 200/Af) 
at the world-tube, for an extraction radius R — 20M. 



Variable 


RateRe 


Rateim 


P 


2.02 


2.01 


J 


2.03 


2.00 


J ,x 


2.04 


1.99 


Q 


2.02 


2.04 


u 


2.02 


2.02 


w 


2.01 


2.04 



We are reporting only first-order convergence rates at future null infinity I + for the Bondi News B and the Weyl 
complonent ^P, but the error is relatively small (0.5% during the late inspiral). From an extraction point of view, these 
errors are smaller than the error in the Cauchy code data and are of little concern. The data are not convergent at 
early time when high frequencies dominate the error. For a thorough analysis of the causes for the first-order accurate 
results and major improvements to the code see [l8j]. 

Figure 2] compares the imaginary and real parts of the (l,m) — (2,2) mode of the Cauchy ip 4 with the complex 
conjugate of the {I, m) = (2, —2) mode of the characteristic "J. We obtain very good amplitude match. Also, we 
observe improved phase agreement as the extraction radius is increased (from 50M to 100M), because the phase error 
in tp4 is reduced with increased extraction radius. 

Figure [5] compares the imaginary and real parts of the (l,m) = (4,4) mode of the Cauchy with the complex 
conjugate of the (I, m) = (4,-4) mode of the characteristic "J. Here we see two effects. First the improved phase 
agreement as Re — ► oo, but also an attenuation of the amplitude due to dissipation of higher-order modes. Also, the 
noise is apparent for the (4, 4) mode. 

Figure [6] compares the amplitudes and the phases between the absolute value of the (l,m) = (2,2) mode of the 
Cauchy ip 4 extracted at R = 50, and the absolute value of the (l,m) — (2, —2) mode of the characteristic "J for the 
same extraction radius, at the highest resolution (N = 200). The difference in amplitude is relatively small, maximum 
0.17% of the Cauchy ip 4 amplitude in the wave zone. 

Figure [7] compares the real part of the (l,m) = (2,-2) modes of the characteristic \& extracted at three different 
extraction radii: R = 20M, R = BOM and R = 100M. The waveform extracted at R = 20 has the biggest amplitude, 
and a very small attenuation of the signal with the radius is observed. 

VI. CONCLUSION 

We have presented here a method for interfacing outer boundary data from a Cauchy evolution with inner boundary 
data for a characteristic evolution so that the waveform can be accurately extracted at infinity. We have demonstrated 
how the PITT null code can be interfaced with the LazEv code, which is a finite-difference BSSN code, to produced 
calibrated waveforms from a binary black-hole inspiral. The extraction interface has been implemented as a thorn in 
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FIG. 4: A plot that compares the phase in the (I = 2, m = 2) mode of the Cauchy t[>4 as calculated using the Null code and 
the Cauchy code. All plots were translated so that the time of the maximum in the amplitude agree. 



0.01 




FIG. 5: A plot that compares the phase in the (2 = 4, m = 4) mode of the Cauchy ip4 as calculated using the Null code and 
the Cauchy code. All plots were translated so that the time of the maximum in the amplitude agree. 



10 




-0.04 



-100 



-50 



0.002 
0.001 


-0.001 



-0.002 
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at 7? = 50, and the (I, m) — (2, —2) mode of the characteristic for the same extraction radius. 
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FIG. 7: A plot of that compares the real part of the (I, m) — (2, —2) mode for the characteristic 9 extracted at three different 
radius (world-tubes). The waveforms are translated such that the maximum of the amplitude corresponds to t/M=0 



the Einstein Computational Toolkit [35[ In this paper we are reporting only preliminary results (see [l8| improvements). 
Although we are aware of the deficiencies in the characteristic waveform extraction tool presented here, there is pressing 
interest from several numerical relativity groups to apply the tool to extract waveforms from binary black-hole inspirals. 
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